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Abstract 



We study a synchronization mechanism, based on one-way coupling of all- 
or-nothing type, applied to coupled map lattices with several different local 
rules. By analyzing the metric and the topological distance between the two 
systems, we found two different regimes: a strong chaos phase in which the 
transition has a directed percolation character and a weak chaos phase in 
which the synchronization transition occurs abruptly. We are able to derive 
some analytical approximations for the location of the transition point and 
the critical properties of the system. We propose to use the characteristics of 
this transition as indicators of the spatial propagation of chaoticity. 
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I. INTRODUCTION 



Recently, the synchronization of chaotic systems has received considerable atten- 
tion milJl,Ei,ll-[Il • Among the many papers on this subject, some of them concern the 
behavior of spatially extended chaotic systems [|3|-|I0| . 



In this paper we study the synchronization properties of two coupled map lattices pi 
when the coupling between them is completely asymmetric, and the synchronization mech- 
anism is of all-or-none kind: denoting one system as the master, and the other as the slave 
one, either an individual map of the slave system is completely synchronized to the corre- 
sponding map of the master system, or it is left free. We study the annealed version of this 
kind of coupling, by choosing at random a fixed fraction of sites to be synchronized at each 
time step, which constitutes the control parameter of the synchronization transition. 

An alternative point of view origins from the problem of the characterization of chaos 
in high dimensional dynamical systems. We assume a rather unusual point of view: instead 
of considering the rate of divergence of the distance between the trajectory of a reference 
system (the master) and a perturbed one (the slave), procedure that leads to the definition 
of the maximum Lyapunov exponent, we measure the "efforts" needed to make the slave 
system coalesce with the master. 

Let us illustrate in detail this approach. The standard technique for studying the chaotic 
properties of such a system is that of measuring the response of the system to a perturbation. 
The conceptual experiment is the following: (1) at a certain time make a copy of the master 
system and slightly vary its state; (2) let the system evolve for a small interval of time and 
(3) measure the distance between the master and the slave; (4) renormalize the distance so 
that it is always small. The logarithmic average of the growth rate of the distance gives 
the maximum Lyapunov exponent. For a small distance, the leading contribution in the 
diverging rate comes from the maximum eigenvalue of the product of the Jacobian of the 
evolution function computed along the trajectory. By considering more than just one vector 
in the tangent space, and keeping them orthogonal, one can obtain the whole spectrum of 
Lyapunov exponents [T^ . 

The Lyapunov spectrum, however, does not describe accurately the process of spatial 
propagation of chaos. We shall illustrate this point by considering diffusively coupled logistic 
map lattices at Ulam point. By analyzing the behavior of the Lyapunov spectrum, one sees 
that for a small coupling the spectrum is almost constant and positive. When increasing the 
coupling between maps, both the positive part of the Lyapunov spectrum and the maximum 
exponent decrease. 

From a different point of view, the chaoticity (and the Lyapunov exponent) of a single 
map can be defined by means of the "efforts" needed to synchronize the slave with the 
master. Let us consider the following "toy" system composed of two simple maps 

= /(*); m 

= (1-P)f{y)+Pf(x). {i> 
where x = x(t) is the master system and y = y(t) is the slave [|14|| ; the prime denotes the 



value of the map at time t + 1. The master system x evolves freely, while the slave y is 
subjected to two opposed contributions: it tends to separate from the master if the map / 
is chaotic, but it is pushed towards x by the parameter p, which represents the "strength" . 
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For a small difference z(t) = x(t) — y(t) one has, 

z' = (l-p)f(x)z (2) 

i.e. the synchronization (z(oo) = 0) occurs for p > p c = 1 — exp(— A), where A is the 
Lyapunov exponent of the (unperturbed) trajectory x(t) of the master system. Thus for 
simple maps the synchronization threshold is related to the chaoticity of the system. We 
would like to extend this concept to spatially extended system, with a diffusive character. 

The basic idea of our approach is to consider Eq. (JTJ) as a mean field description of a 
stochastic process which, with probability p, sets each individual map of the slave system 
to the value of the corresponding map of the master one ( "pinching" ) . We assume that the 
synchronized state is absorbing, i.e. once a patch of the slave system is in the same state 
of the master one, desynchronization can occour only at borders, without bubbling. The 
completely synchronized state, even if unstable, cannot be exited. A discussion about the 
robustness of this synchronization mechansm is deferred to the last Section. 

This method can be applied also to more exotic dynamical systems, such as cellular 
automata ]T3| . For such systems the usual chaotic indicators cannot be easily computed, 
but they exhibit spatial propagation of disorder. 

It is clear that for p = 1 the system synchronizes in just one time step. However, as 
we shall show in details in Section 3, a critical value p = p c < 1 does exist, above which 
the system synchronizes regardless of its chaoticity. This threshold is related to a directed 
percolation phase transition. 

On the other hand, some systems can synchronize for p = p* < p c , because the propaga- 
tion rate of the differences between the master and the slave is reduced by the inactivation 
of the degrees of freedom due to pinching. As a trivial illustration of this process, a lattice 
of uncoupled maps synchronizes in the long time limit for all strengths p > p* = 0. This syn- 
chronization threshold is an indicator of the spatial propagation of chaoticity in the system. 
Other informations come from the dynamics of the transition. 

The sketch of this paper is the following: in the next Section we describe precisely the 
model we use and introduce the observables, in Section 3 the synchronization transition is 
analyzed, and the connections with the directed percolation problem discussed. In Section 4 
we present the phase diagrams for the synchronization transition for several well-known 
maps, and finally in the last Section we discuss our main conclusions, including possible 
extensions and applications of the synchronization mechanism, and the relationships with 
other models found in literature. 

II. THE STOCHASTIC SYNCHRONIZATION MECHANISM 

The state of the master system at a given time t is denoted as X 1 e [0, 1]^, and a 
component of X (a single map) is indicated as X\i = X{, i = 1, . . . , N. The dynamics is 
defined as 

= gM-iAA+i) = (1 - 2e)/(xJ) + e(/(a£.i) + f{<-i)\ (3) 

with periodic boundary conditions and < e < 1/2. It can be considered as the discretiza- 
tion (in time and space) of a reaction-diffusion system, the e-coupling coming from the 
Laplacian operator. The dynamics of the slave system Y is given by 
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Vi = (i - nip))M-i, yl Vl+i) + rMgeixU,^, xl +1 ), (4) 

where r\{p) is a random variable which assumes the value one with probability p, and zero 
otherwise. In other words, at each time step a fraction p of maps in the slave system is set 
to the same value of the corresponding map in the master system. 
In vectorial notation one can write for the master system 

X t+l = G £ (X') = (/ + eA)F(X t ), 
where A is the discrete Laplacian 

AX\ij = (Sij-i + 6ij+i — 2Sij)xj 
and F(X) is a diagonal operator F(X)\ij = f(xj)Sij. For the slave system one has 

Y t+1 = s\p)G £ (Y t ) + S\p)G £ (X t ) 

where S(p) is a random diagonal matrix having a fraction p of diagonal elements equal to 
one and all others equal to zero, S(pY\ij = r t j(p)S i j, and S(p) = I — S(p). 

We introduce also the difference system Z = X — Y, whose evolution rule is 

Z t+1 = s\p)(G £ (X t )-G £ (X t -Z t )) 

If the difference field is uniformly small, i.e. for max(|4l) — ► 0, one has 

Z t+1 = s\p)J £ (X t )Z\ 

where J £ (X) is the Jacobian of the evolution function 

J c (X') = (I + sA)^l 

and 



dF(X) 



dX 



df(x) 
= 5 -^x~ 



For all kinds of numerical computations there is a limit to the precision below which 
two numbers become indistinguishable. Since this limit depends on the magnitude of the 
numbers, it is very hard to control its effects. In this perspective, we introduce a threshold 
t on the precision, by imposing that if \xi — Ui\ < r then yi = X{. In this way we can 
study the sensitivity of the results on r, and eventually perform the limit r — > 0. We have 
checked that our asymptotic results are independent of r, at least for r smaller than 10~ 6 . In 
the following we shall neglect to indicate the truncation operation for the ease of notation, 
except when explicitly needed. Due to the precision threshold r, the difference field z\ is set 
to zero if z\ < r. 

An alternative way of computing the evolution of the system, which will be useful in the 
following, is given by the following procedure: 
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1. Consider three stacked two-dimensional lattices of size N x T, and label one direction 
as space i, i G [l,iV] and the other one as time t, t G [0, T\. This ensemble can be 
considered composed by three layers: one that contains the Boolean numbers rj, and 
two containing the real numbers x\ and y\. 

2. Fill up the r\ layer with zeros and ones so that the probability of having r\ = 1 is p; 
this layer will be named the quenched field. 

3. Fill up the x® and y® rows at random; they will be the initial conditions for the X and 
Y lattice maps. 

4. Iterate the applications (|3|) and to fill up the X and Y layers. 

In this way we define an out of equilibrium statistical system, with p as a control pa- 
rameter. It is assumed the limit N —>■ oo and T — > oo and the average over the quenched 
field. The degree of synchronization of the system at a given time t can be measured by the 
(metric) distance 

1 N-l 

C l = — V \z l \ 

iv i=0 

whose asymptotic value for a given probability p will be denoted as CO )- 
We introduce also the field M as 



and the topological distance 



if z\ = 0; 

1 otherwise, 



1 N 
P l = -Th\ 
F N j-t 1 



which measures the fraction of non-synchronized sites in the system. The asymptotic value 
of the topological distance will be denoted as p(p). 

We shall study the synchronization transition for the following maps f(x) of the unit 
interval x G [0, 1]: 

1. The generalized Bernoulli shift f(x) = [aa;]modl with a slope a greater than one. The 
Lyapunov exponent of the single map is simply ln(a); this is also the value of the 
maximum Lyapunov exponent for a diffusively coupled lattice, regardless of e. 

2. The quenched random map that assumes a different (random) value for each different x. 
It can be considered equivalent to the Bernoulli shift in the limit a — * oo, and thus the 
map with the highest degree of chaoticity. This map is everywhere non-differentiable. 

3. The logistic map f(x) = ax(l — x), which is chaotic for 3.57 . . . < a < 4. 

4. The generalized tent map 



with a > 1. 
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III. THE SYNCHRONIZATION TRANSITION 



If the maps are uncoupled (e = 0), any value of p greater than zero is sufficient to 
synchronize the system in the long time limit, regardless of the chaoticity of the single map. 
For coupled systems, however, the chaoticity of the map (temporal chaos) contributes to the 
spatial chaoticity of the system. 

For illustration purposes, we show in Fig. |l] the behavior of the metric distance ((p) (left) 
and the topological distance p(p) (right) for Bernoulli maps with e = 1/3 and different values 
of a. The topological distance exhibits a sharp transition for a < a c ~ 2, and a smooth 
transition for a > a c ; all curves superimposes to a universal curve far from the transition 
point. The metric distance always exhibits a smooth transition. This scenario is generic for 
all kinds of maps and couplings, except that some maps (notably the logistic map) never 
exhibit the smooth transition of the topological distance in the allowed range of values of 
the parameter a. 

For the random map it is quite easy to understand the origin of the universal curve: 
since even a small distance is amplified in one time step to a random value, the difference 
z\ is greater than if z 1 ^ 1 > on some of the neighbors (\i — j\ < 1) at the previous time 
step, and the quenched field r\ is equal to 0. 

This defines a directed site percolation problem [15] (formation of a percolating cluster 



of "wet" sites along the time direction), with control parameter 1 — p. A site of coordinate 
(i,t) is said to be wet if rj(p) = and it is connected at least to one neighboring wet site 
r* _1 (p) = (\i — j\ < 1) site at time t — 1. All sites in the first (t = 0) row are supposed 
to be connected to an external wet site. We shall denote the ensemble of wet sites with the 
name wet cluster; an instance of a wet cluster for p = 0.45 is shown in Fig. [|. Notice that 
the usual directed percolation probability p( DP ) is given here by p( DP ) = 1 — p. 

If j/ DP ) < pp p ) the wet cluster does not percolate in the time direction, and the system 
will synchronize despite any chaoticity of the single map. The r\ = sites disconnected from 
the percolation cluster do not influence the synchronization, since the synchronized state is 
locally absorbing. Thus, for the random map and e > 0, the invariant curve is simply the 
curve of the asymptotic density m(p) of wet sites for the directed percolation problem. In 
the vicinity of the transition, m(p) behaves as 

™>(p) ~ (p c - P) P , 

where (3 ~ 0.26(1) is the magnetic critical exponents of the 1 + 1 dimensional directed per- 
colation problem, and p c — 1 — p^P^ — 0.460(2) for this lattice with connectivity three [16 



In the generic case, one has to study what happens on the wet cluster. If the evolution 
of the system is expanding in the difference space, and the wet cluster percolates in the 
time direction, then the asymptotic topological distance is greater than zero and vanishes 
when the density of wet sites does. The expansion rate in the difference space depends on 
the average number of connections between wet sites in the cluster, an exception being the 
random map, which expands to a random value regardless of the initial difference. 

On the other hand, it may happen that for some value of p (for a given percolating wet 
cluster) the dynamics in the difference space is contracting, so that the difference field z\ 
will eventually become less than r and thus set to zero. As an illustration, let us consider 
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the Bernoulli shift for e = 1/3 and only one connection in average. The problem reduces to 
the computation of the expansion rate in a one dimensional chain of asymmetrically coupled 
maps, and the threshold for the synchronization is a = 3. 

This eventual contraction mechanism implies that the cluster of the sites for which h\ = 1 
(the difference cluster) is generally strictly included into the percolation cluster. 

We classify the case in which the difference cluster always survives when the percolation 
cluster spans the lattice with the name strong spatiotemporal chaos, otherwise we have the 
weak spatiotemporal chaos. We shall denote with a star the values of the parameters for 
which the synchronization transition occurs, distinguishing the weak and strong cases by 
the index w or s, respectively. 

From the results of our numerical simulations, we think that the following Ansatz for 
the topological distance p(p) holds: the strong chaos transition always belong to the di- 
rected percolation universality class, while the weak chaos transition is always of first order 
character. This statement is equivalent to the assumption that the fraction of sites in the 
percolation cluster that do not belong to the difference cluster is either vanishing (strong 
chaos) or order one (weak chaos) in the thermodynamic limit. If noise is allowed to desyn- 
chronize the system, the directed percolation character is lost. 

This Ansatz can be illustrated by plotting the ratio between the topological distance 
p(p) and the density of wet sites m(p) as shown in Fig. ^ for the logistic map and e = 1/3. 
One can see that this ratio maintains almost constant up to the transition point. We have 
checked that this behavior hold for several values of the slope a and the coupling e for the 
three maps. 

It is interesting to examine the behavior of the difference z\ near the transition point. 
In the strong chaos phase, the value of the difference z\ for non-synchronized sites is large, 
since the synchronization if forced by the percolation mechanism. On the other hand, in 
the weak chaos phase, the value of the difference z\ for non-synchronized sites is small in 
average, even though localized bubbling of large difference can be observed [17 ]. 

In the following we study numerically the phase boundary between the strong and weak 
chaos for the Bernoulli shift and the logistic and tent map, and we obtain some analytical 
approximation for it. 



IV. PHASE DIAGRAMS 

The boundary between the weakly and strongly chaotic phases can be defined considering 
the synchronization mechanism on the critical percolation cluster, i.e. the cluster of wet sites 
at p* = minp* = maxp^ = 1— pp p ). As illustrated in Appendix A, the critical wet cluster is 
included in all percolating wet clusters, so it can be considered to be the smallest percolating 
cluster of wet sites, and thus the one most unfavorable to the proagation of the difference, 
still presenting spanning paths over which chaos can propagate. The algorithm for the 
generation of the critical cluster is also described in Appendix A. 

The phase diagrams for the Bernoulli shift and the tent map are shown in Fig. [3]; the 
logistic map never exhibits the strongly chaotic phase. The case e — 1/2 is special, since the 
connectivity changes, and in this case the lattice corresponds to that of the Domany-Kinzel 
model |I8|, for which 1 — p^" 1 = 0.3. 
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Since this boundary phase is defined on a critical percolation cluster, one can assume that 
the average expansion rate, given by the sum of all paths in the wet cluster, is dominated 
by a single path. Let us assume that this path contains x vertical steps (of weight 1 — 2e) 
and y oblique steps (of weight e). For the Bernoulli shift, all steps carry an expansion rate 
a and thus one obtains for the Lyapunov multiplier p c at the phase boundary 

p c ~ a x+y (l - 2e) x e y ~ 1 

and thus 

const 

a — ~, ; i 

(1 - 2e) Q i£ Q 2 

with 

oi\ + a 2 = 1. (5) 

Indeed, the numerical data presented in Fig. ^| does support a correspondence of this form. 
However, they are too rough to allow the precise computation of the exponents, which 
anyhow do not seem to always obey to Eq. (|5|). 

We computed numerically the transition point p* and the critical exponents (3 and u± 
for the topological distance p(p) for the Bernoulli shift and the tent map by means of the 
scaling relation 

~) ■ p(a^(p-p*)+p*,— ) (6) 

that holds in the thermodynamic limit iV — > oo. Here a is an arbitrary time scaling factor. 
We found p c = 0.460(2), (3 = 0.26(1) and u± = 1.75(5). The results are consistent with the 
hypothesis that the phenomenon belongs to the directed percolation universality class, as 
expected. 

The strong chaos transition is characterized by the expanding dynamics on the perco- 
lation cluster. Thus, in general, at the transition point the value of a non-zero difference 
z\ is not small. On the contrary, the weak chaos transition is characterized by the average 
vanishing of z\. Up to the boundary between the two regimes the distance field Zi is still 
small and a linear approximation can be used. The synchronization mechanism is related to 
the average number of links between wet sites, k(p). Since the wetting is forcing (once that 
a site has been wet, it cannot be "unwet"), one can simply consider the process for a single 
site, obtaining n(p) = 3(1 —p). 

By imposing that the distance should be marginally expanding, and assuming that the 
difference field Zi is almost constant, we get for the Bernoulli shift and e = 1/3, 

1 

3<«(Pw) = <C0--Pw) = 1 

and thus a* w = 1/(1 -p* w ). 

For a generic map, with positive and negative slopes, cancellation effects can be present. 
We need an indicator about the constancy of the sign of the derivative over the interaction 
range of the diffusion operator. We consider the quantity 
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f'( Xi ) + f(x i+1 ) 
Vi- 2 

which is analogous to the local Lyapunov multiplier, and 

/N-l \Vf 

v = y II 

as a global indicator. For the Bernoulli map this indicator corresponds to the slope a. 

Thus, for e = 1/3, the condition for observing the weak chaos synchronization transition 
is given in this approximation by 

3 • »7« ■ ) = rf w {l - p* w ) = 1. 

We report in Fig. [5| the behavior of ^(1 — p^) as a function of the slope a for Bernoulli 
shift and tent maps, together with the mean field approximation. The discrepancies from 
this approximation originate from the assumption of a uniformly vanishing difference field 
Zi, which is not fulfilled even at the transition point. 

Let us turn now to the metric distance (. For random maps even an infinitesimally small 
distance is amplified in one step to a random value. Thus, we have for the metric distance 

C(*) = (WW), 

and ( has the same critical behavior of p, denoting a certain degree of universality in this 
synchronization transition. 

We have studied the behavior of the ratio C(p)/p(p) as a function of p for various values 
of the slope a for the Bernoulli shift and the tent map, and e — 1/3. The results, reported in 
Fig. ^ show that for the Bernoulli shift the metric distance is independent of p and of slope 
a far from the synchronization transition. Deviations from this behavior near the transition 
point vanish as the slope a become greater than a*. 

For tent map the ratio ((p) / p(p) is independent of p far for the transition but it depends 
on the slope a which determines the distribution of Zj. Again, for a ^> a*, the metric distance 
is almost constant from p also near the transition point. 

The case e ^ 1/3 is more difficult to analyze, since the asymmetric couplings cannot 
be easily mapped onto a statistical problem. However, for very small e, the expansion is 
dominated by the exponential growth (with average rate A) along the vertical link of the 
largest difference. A very rough mean field description could be the following. Let us assume 
that at a certain time there is essentially only one site i with a non null difference Z{. This 
difference grows exponentially at rate A, and propagates to the neighboring sites at rate 
e. After an average time 1/p the difference at site i is set to zero by the synchronization 
mechanism, and after time 2/p only one of the neighboring has a nonzero difference, thus 
closing the cycle. Thus, the average expansion rate at the synchronization threshold p c is 

e exp(2A/p c ) ~ const, (7) 

where A is the Lyapunov exponent of the uncoupled map, and we have neglected non- 
exponential prefactors. In Fig. [7] we show the results of one simulation for the three maps 
studied in the article, with the parameters chosen so to have A = ln(2), and for Bernoulli 
maps with various slopes. One can see that Eq. (^) is verified for small e, except finite size 
and time effects. 
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V. CONCLUSIONS 



We have studied the synchronization transition between two chains of diffusively coupled 
chaotic maps, induced by the inactivation of degrees of freedom in the difference space with 
a probability p. We have found that two different regimes can be defined: the strong chaos 
regime for which the dynamics of the transition is dominated by the directed percolation 
transition, and a weak chaos regime in which the system synchronizes in presence of spanning 
paths along which the difference could in principle survive. We have been able to present 
some analytical approximation of the transition point and its critical properties. 

The character of the transition and the critical value of the parameter p are proposed as 
indicators of spatial propagation of chaoticity, which can complement the usual Lyapunov 
description. These indicators does not rely on the existence of a tangent space or exponential 
growth, so they can be applied to a broader class of systems, like non differentiable maps or 
cellular automata JT3|], and system presenting stable chaoscitestablechaos. 



Our approach can be considered as the annealed version of models that exhibit a synchro- 
nization transition, presented in some recently appeared papers. First of all, let us consider 
the synchronization mechanism proposed by Pecora and Carroll [|1]. In their numerical and 
experimental set-up they studied the behavior of the distance between two chaotic oscilla- 
tors, when part of the degrees of freedom of one of them is set equal to the corresponding 
degrees of freedom of the other. Extending this mechanism to spatially extended systems, 
one has a coupling similar to ours, but with quenched disorder (the coupling degrees of 
freedom). We did not perform the study of quenched disorder in diffusively coupled maps, 
since it implies longer spatial couplings in order to avoid the formation of walls. 

Another similar system was studied by Fahy and Hammann [^C]. Their subject was 
an ensemble of non-interacting particles in a chaotic potential. At fixed time intervals the 
velocities of the particles were all set equals to a Gaussian sample. If the free-fly time was 
small enough, all trajectories collapse into one. We performed preliminaries simulations (not 
reported here) on a modified system, in which one reference particle followed an unperturbed 
trajectory, while a replica had its velocity set equal to that of the reference at time intervals 
r. Indeed, we observed a synchronization transition for small enough r. By observing the 
system stroboscopically at intervals r we can substitute the continuous dynamics with a 
map. In this case the time r controls the chaoticity of the map, while the synchronization 
mechanism is similar to that of Pecora and Carroll and thus to a quenched disorder. Fahy 
and Hammann also checked that the synchronization transition occurs if the position of the 
particles is set equal, instead of the velocity. 

The synchronization mechanism studied in this article is quite particular, since it implies 
a complete collapse of the distance between the master and the slave. The strong chaos 
synchronized phase is not stable with respect to the inclusion of desynchronizing effects, 
such as noise or non perfect identity of parameters in the master and slave systems. On the 
other hand, the weak chaos transition is ruled by the exponential shrinking of difference, i.e. 
by negative Lyapunov exponents for the difference. Thus, it is expected that this transition 
is robust with respect to weak desynchronization effects, and not to depend on the threshold 
r. The modification of the observed phenomena in presence of noise will be the subject of 
a future work. In the present version, the synchronization transition can be considered 
as a mathematical tool for the definition of quantities related to the spatial propagation 
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of chaoticity. The most natural systems to which this method can be applied are those 
presenting stable chaos (jnj, i.e. irregular behavior in the presence of negative Lyapunov 
spectra. For these systems an eventual small noise if wiped out by the contracting dynamics 
on small scales. Similar noise-free systems are those that can be approximated by cellular 
automata models 113 



ACKNOWLEDGEMENTS 

We thank S. Ruffo, R. Rechtman and M. Bezzi for fruitful discussions. Part of the sim- 
ulations have been performed on the CRAY-T3E of the CINECA center using the INFM 
facilities (Progetto Calcolo Parallelo). L.B. and P.P. thank the Dipartimento di Matemat- 
ica Applicata for kind hospitality. Part of this work was performed during the workshop 
"Complexity and Chaos" at ISI Foundation, Torino, Italy. 



APPENDIX A 

The generation of the critical cluster can be done in a very efficient way, using a modifi- 
cation of the invasion percolation algorithm |2T . 



The random bit rj(p) is obtained by comparing a random number Rj, whose distribution 
probability is constant in the unit interval, with p. If Rj < p then rj — 1, else rj = 0. 

Let us consider a lattice with the same geometry of the percolation one and assign to 
each site a random number Rj. The idea now is that of lowering p (starting from p = 1) 
until the cluster of rj = sites spans the lattice. A nice property of the directed site 
percolation problem (which is related to the "forcing" character of wetting) is that, given 
the random number Rj, the percolation cluster for p\ is included in the percolation cluster 
for p2 if pi < P2- Thus, one has simply to choose the site with the highest Rj on the border 
of percolation cluster (i.e. among the sites with rj = l connected to some r* _1 = site, with 
\i — j\ = 1). This maximum value will become the new estimate of p, and the percolation 
cluster is enhanced to include all sites with Rj > p. 

This procedure can be easily performed by keeping the values of the sites in the border 
in an ordered linked list, assuming that all sites in the t = row are connected to a wet site. 
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FIG. 2. An example of a wet cluster for p = 0.45. 
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FIG. 3. Percolative phase diagram for Bernoulli (a) and tent (b) maps. The line represent a 
power law fitting, as described in the text. Average over 5 runs, N = 100, T = 4000. 
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FIG. 4. The ratio p(p)/m(p) vs. p for a chain of 1000 logistic maps, T = 3000, slope a = 4. 
and £ = 1/3. 
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FIG. 5. The relation between the slope a and the weak chaos synchronization threshold 
77^/(1 — p^) for a chain of Bernoulli shift (a) and tent maps (b). Average over 4 runs, N = 400, 
T = 8000, e = 1/3. 
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FIG. 6. The dependence of the ratio between the metric distance Q{p) and the topological 
distance p(p) as a function of p, for the Bernoulli shift (a) and tent map (b) and different slopes a. 
Average over 4 runs, N = 100, T = 1000. 
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FIG. 7. Relation between a small coupling e and the synchronization threshold p^. Part (a): 
the three sets correspond to the Bernoulli shift with slope a = 2 (stars), tent map with slope a = 2 
(diamonds) and logistic maps with a = 4 (squares). For all three maps the maximum Lyapunov 
exponent for e = is A = m(2); the dashed line correspods to the law = ln(e)/(21n(2)). Part 
(b): Bernoulli shift for a = 2 (circle), a = 3 (triangles) and a = 4 (squares); the dashed lines have 
slope 1/(2 ln(a)). Data from one simulation with N = 500 and T = 2000. 
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